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We study the semileptonic decays of Df, D + , and D° mesons into the light scalar mesons [/o(500), 

Ao(800), /o(980), and ao(980)] and the light vector mesons [p(770), w(782), JT*(892), and 0(1020)]. 

With the help of a chiral unitarity approach in coupled channels, we compute the branching fractions 
for scalar meson processes of the semileptonic D decays in a simple way. Using current known values 
of the branching fractions, we make predictions for the branching fractions of the semileptonic decay 
modes with other scalar and vector mesons. Furthermore, we calculate the 7r + 7r _ , ng, -kK, and 
K + K~ invariant mass distributions in the semileptonic decays of D mesons, which will help us 


clarify the nature of the light scalar mesons. 
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I. INTRODUCTION 

The recent experimental situation in hadron physics 
enables us to utilize huge amounts of data on heavy 
hadrons, which contain charm or bottom quark(s), for 
the investigation of hadron structures. Especially, de¬ 
cay properties of heavy mesons can shed more light on 
the nature of the light scalar mesons [/o(500), A'q(800), 
/o(980), and ao(980)], which has been a hot topic in 
hadron physics [1]. For instance, the decay —»■ 

J/ipTT + n~ has been experimentally measured in Refs. [2- 
6] for the study of the /o(500) and /o(980) resonances, 
and they observed a pronounced peak for the /o(980) 
while no evident signal was found for the /o(500). Then 
a theoretical study [7] followed the experiments and re¬ 
produced ratios of experimental branching fractions at a 
quantitative level, pointing out that J /ip + (ss) produc¬ 
tion in the B° s decay and a hadronization of ss to I\K are 
essential to understand the branching fractions of the B° 
decay into J/0’/o(98O). In the theoretical study, the final 
state interaction between two pseudoscalar mesons is cal¬ 
culated with the so-called chiral unitary approach [8-16], 
in which the light scalar mesons are obtained as dynam¬ 
ically generated resonances, and it is concluded that the 
/o(980) has a substantial fraction of the strange quarks. 
The same hadronization scheme has been employed in 
theoretical studies in Refs. [7, 17-19]. 

In this paper, we consider the semileptonic decay of 
D -a hadron(s) + l + vi , extending a discussion for the 
semileptonic B decays into D* 0 (2317) and Dq(2400) reso¬ 
nances in Ref. [19]. The semileptonic D decays have been 
experimentally investigated in, e.g., BES [20, 21], FO¬ 
CUS [22, 23], BaBar [24, 25], and CLEO [26-30], Here, 
in order to grasp how the semileptonic decay takes place, 
let us consider the Df meson. Since the constituent 
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quark component of Df is cs, we expect a Cabibbo fa¬ 
vored semileptonic decay of c —)• sl + vi and hence the 
decay Df —> (ss) l + vi with ss being the vector meson 
0(1020), which is depicted in Fig. 1(a). Actually this 
semileptonic decay mode has been observed in experi¬ 
ments, and its branching fraction to the total decay width 
is B[Df -a 0(1020) e+ u e ] = 2.49±0.14% [1] (see Table I, 
in which we list branching fractions for the semileptonic 
decays of Df, D + , and D° reported by the Particle Data 
Group). On the other hand, we cannot straightforwardly 
extend the discussion to the scalar meson productions 
in the final state of the semileptonic decays, since the 
structure of the scalar mesons, whether qq or some ex¬ 
otic one, is still controversial. In this study we consider 
the production of the /o(980) or / 0 (500) as dynamically 
generated resonances in the semileptonic Df decay, so we 
have to introduce an extra qq pair to make a hadroniza- 




FIG. 1: (a) Semileptonic decay of D+ into l + isi and a primary 
ss pair, (b) Semileptonic decay of Df into l + vi and two 
pseudoscalar mesons P with a hadronization. 
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TABLE I: Branching fractions for the semileptonic decays of 
Df, D + , and D° reported by the Particle Data Group [1]. In 
this Table we only show decay inodes relevant to this study. 


D+ 

Mean life [s] 


(500 ± 7) x 10" 15 

1020)e + i/ e ] 


(2.49 ±0.14) x 10 -2 

B[uj(782)e + v e ] 


< 2.0 x 10“ 3 

B[K* (892)° e + v e ] 


(1.8 ±0.7) x 10~ 3 

6[/o(980)e+i/ e , /o(980) —I w 

■+7T ] 

(2.00 ± 0.32) x 10' 3 


D+ 


Mean life [s] 


(1040 ±7) x 10” 15 

B[K*(892)°e + i/ e , K* (892)° 

-s- K~tt+] 

(3.68 ±0.10) x 10” 2 

B[(K~n + ) s - wave e + iy e ] 


(2.32 ±0.10) x 10” 3 

B[K*{ 892) ( V+iv, 7?*(892)° 

-s- K~tt + ] 

(3.52 ±0.10) x 10” 2 

B[p{770)°e + v e } 


(2.18±g;|5) x 10~ 3 

B[p( 770)V^] 


(2.4 ±0.4) x 10^ 3 

B[uj(782)e + Ve] 


(1.82 ±0.19) x 10” 3 

B[<j>{ 1020)e + i/ e ] 


< 9 x 10“ s 


D° 


Mean life [s] 


(410.1 ± 1.5) x 10- 15 

B[K* (892) ~ e + v e \ 


(2.16 ±0.16) x 10~ 2 

B[I<*{ 892)-/i+i^] 


(1.90 ±0.24) x 10" 2 

B[A' - 7T°e + i'e] 


(1-6 tH) x 10- 2 

Z3[A'°7T _ e + ^e] 


(2.7±o. 7) x 10" 2 

B[p{770)-e+u e ] 


(1.9 ±0.4) x 10~ 3 


tion as shown in Fig. 1(b). The introduction of an ex¬ 
tra qq pair to make a hadronization has been performed 
in Refs. [7, 17-19]. In this study we apply the same 
method of the hadronization to the semileptonic decays 
of D mesons so as to investigate the nature of the light 
scalar mesons. 

Utilizing the semileptonic decay of a heavy hadron pro¬ 
vides us with two advantages when we investigate the 
internal structure of hadrons in the final state of the 
semileptonic decay. First, Cabibbo favored and sup¬ 
pressed processes enable us to specify flavors of quarks 
contained in final state hadrons. Second, the semilep¬ 
tonic decay of the heavy hadron to two light hadrons 
+l + vi brings a suitable condition to measure effects of 
the final state interaction of the two light hadrons, since 
the leptons and hadrons in the final state interact with 
each other only weakly. 

Theoretical work on the issues of the semileptonic D 
decays is already available. In Ref. [31], using QCD sum 
rules, the Df and D + semileptonic decays into /o(980) 
are considered concluding that the importance of up and 
down quarks in the / 0 ( 980) is not negligible. In Ref. [32] 
the D+ —>• /o(980)e + r' e reaction is analyzed from the 


point of view of the /o(980) being a qq state, concluding 
that ss component of the /o(980) may not be dominant. 
In Ref. [33] the Df —> n + -r:~e + u e reaction is studied con¬ 
cluding that it supports the dominant four quark nature 
of the /o(500) and / 0 ( 980). Similar conclusions about 
the four quark nature of the scalar mesons are reached 
in the work of [34, 35]. Research along the same line is 
done in Ref. [36], looking for likely reasonable ratios that 
would help distinguish between the two and four quark 
structure of the scalar mesons. 

Another line of research is done using light-cone sum 
rules to evaluate the form factors appearing in the pro¬ 
cess [37]. This line of research is applied in many related 
processes, rare decays like B s tt + tt~1 + 1~ in [38], B s —> 
in [39], B ^ —> J/i/jTr + n^ and B s 7r + 7r _ /i + /r~ 
decays in [40], or semileptonic decays [41-43]. In some 
cases the meson final state interaction is further imple¬ 
mented using the Omnes representation [37, 42], while in 
other cases Breit-Wigner or Flatte structures are imple¬ 
mented and parametrized to account for the resonances 
observed in the experiment. 

In contrast to these pictures, in the present study we 
treat the scalar mesons as dynamically generated reso¬ 
nances from two pseudoscalar mesons in the so-called chi¬ 
ral unitary approach. Then we describe the semileptonic 
decays of D mesons in an economical way for hadroniza¬ 
tion as done in Refs. [7, 17 19]. 

This paper is organized as follows. In Sec. II we for¬ 
mulate the semileptonic decay widths of Df, D + , and 
D° into the light scalar and vector mesons and give our 
model of the hadronization. We also calculate meson- 
meson scattering amplitudes to generate dynamically the 
scalar mesons. In Sec. Ill we show our numerical results 
of the semileptonic decay widths of £>+, D + , and D°. 
We predict branching fractions which are not reported 
by the Particle Data Group and show invariant mass dis¬ 
tributions of the two pseudoscalar mesons from the scalar 
and vector mesons. Section IV is devoted to drawing the 
conclusion of this study. 


II. FORMULATION 

In this section we formulate the semileptonic decay 
widths of Df, D + , and D° into light scalar and vector 
mesons: 


Df, D+, D° -> 


Sl+v u 

Vl+vi, 


S -A PP, 


(1) 


where S, V, and P represent the light scalar, vector, and 
pseudoscalar mesons, respectively, and the lepton flavor 
l can be e and /i. Explicit decay modes are listed in Ta¬ 
ble II. In order to formulate the decay width, we consider 
first the semileptonic decay amplitudes and widths in 
Section IIA and next hadronizations into scalar and vec¬ 
tor mesons in Section IIB. Scattering amplitudes of two 
pseudoscalar mesons are then constructed in the chiral 
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TABLE II: Semileptonic decay modes of Df, D + , and D° con¬ 
sidered in this study. The lepton flavor l is e and p. We also 
specify Cabibbo favored/suppressed process for each decay 
mode; the semileptonic decay into two pseudoscalar mesons 
is judged with the discussions given in Sec. II B. 


D+ 

0(1020) l + vi 


favored 

K* (892)° Z + vi 


suppressed 

7T + 7r _ / + Vl 


favored 

K+K~ 1+ ui 


favored 

7T“ K + l + Vl 

D+ 

suppressed 

A*(892)° l + vi 


favored 

p( 770)° 1+ pi 


suppressed 

oj( 782) 1+ pi 


suppressed 

+ 

1 

' + 


suppressed 

7r°?7 1 + Pi 


suppressed 

K+K~ 1+ pi 


suppressed 

n + K~ l + pi 

D° 

favored 

I<* (892)" 1+ pi 


favored 

O 

1 

* + 


suppressed 

7T p l + Pi 


suppressed 

K°K~ 1+ pi 


suppressed 

7 t~K°1 + pi 


favored 


unitary approach for the description of the scalar mesons 
in Section IIC. Throughout this study we assume isospin 
symmetry for light hadrons. 


A. Amplitudes and widths of semileptonic D 
decays 

In general, we can express the decay amplitude of D —> 
hadron(s) + l + vi, Tb, by using the propagator of the W 
boson and its couplings to leptons and quarks, which can 
be replaced with the Fermi coupling constant Gp. At this 
stage we do not fix the number of the final state hadrons. 
In a similar manner to the formulation in Ref. [19], the 
explicit form of Tjj becomes 

Td = -i^%L a Q a X Vhad- (2) 

v 2 

The factor Vhad consists of the wave function of quarks in¬ 
side the D meson, the hadronization contribution in the 
final state, and the Cabibbo-Kobayashi -Maskawa ma¬ 
trix element for the transition from the charm to a light 
quark. The explicit form of Vhad will be determined in 
the next subsection. The lepton and quark parts of the 


W boson couplings are defined as. 

L = Ujj'y (1 75)^L Qa — 75)^c? (3) 

respectively, where w„, Vi, u q , and u c are the Dirac 
spinors corresponding to the neutrino, lepton l + , light 
quark q , and charm quark, respectively. 

Let us now calculate the squared amplitude for the 
semileptonic D decay widths, in which we average (sum) 
the polarizations of the initial-state quarks (final state 
leptons and quarks). Therefore, in terms of the ampli¬ 
tude in Eq. (2), we can obtain the squared decay ampli¬ 
tude as 

I ]T |r D | 2 = !2dhT iced 2 w 

pol pol 


where the factor 1/2 comes from the average of the charm 
quark polarization in the initial state. We can further 
calculate the lepton and quark parts in the amplitude (3), 
by using the conventions of the Dirac spinors and traces 
of Dirac 7 matrices summarized in Appendix A, which 
lead to 

L a L t/3 =tr 

pol 

_ 2 pfpt + pZPi - Pi ■ Vv9 afi + ie a0pcr PipPv<7 
mim v 

( 5 ) 




■75 


)7 /3 - 


m v 


2m v 


where pi and p u (mi and m v ) are momenta (masses) of 
the lepton l + and neutrino, respectively, and 


^2 QaQl = tr 



75) 


j>c + m e 
2 TO C 


(1 + 75)7/3 


i>q + m q 
2 m q 


2PcaPqp + PqaPcp ~ Pc ' Pgffa/3 + i^appaP P c Pq 
m c m q 

( 6 ) 


with the momenta (masses) of the charm and light 
quarks, p c and p q (m c and m q ), respectively . 1 Then with 
a straightforward calculation we have 


Ei iQ ^i 2 


16 (pi -Pc)(Pu - Pq) 
mim v m c m q 


( 7 ) 


Now let us rewrite the momenta of quarks by using those 
of hadrons in the following manner: 

P^_ = P^^ (8) 

m c mo ’ m q thr 


1 The momentum p q is for a quark in the primary qq pair after 
the W boson emission, which means that the momentum p q is 
carried by the constituent quark. Accordingly, m q is the mass of 
the constituent quark rather than of the current quark. In this 
sense, m q respects the flavor SU(3) symmetry. 
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where we have neglected the relative internal momenta 
of the quarks, which are typically small compared to the 
masses of quarks. Here uid and m R (pd and p R ) are 
the masses (momenta) of the D and R = S, V mesons, 
respectively. With these translations the square of L a Q a 
with polarization summation becomes 


put it out of the integral. In this study, however, W a d 
will be simply constructed, so that this will not depend 
on nor the angle, as we will see in the next sub¬ 

section. Furthermore, the integral of the solid angle 
is performed in the Iv rest frame as [19] 


yu.Q i2 _ 16 {pi ■ pp){p v ■ Pr) 

, “ mim v m D m R 

pol 

Therefore, we obtain the squared decay amplitude as: 
4|G F I4 ad | 2 


(9) 


Je i^i 2 = 


pol 


mim u rriDmR 


■{pi ■ Pd){Pv ■ Pr)- (10) 


J dVL v [pi ■ p D ){Pv -Pv) 

= J dQ. v {EiE D + p v -p D ){E v E v -p u -p D ) 

~ ~ ~ ~ 47T 

= AttEiE u E d E v - —\Pv\~\Pd \ 2 

= *[M ( £>] 2 (e d E v - l -\p D \ 2 X 


(15) 


With the above squared amplitude we can compute 
the decay width. We will be interested in two types 
of decays: three-body decays for vector mesons such as 
Df 0(1020) e + v e , and four-body decays for scalar 
mesons constructed from two pseudoscalar mesons such 
as ot —> 7T + 7r e + v e . As it will be seen, both decay 
types can be described by the amplitude Td with differ¬ 
ent assumptions for Wad: V^ad an< ^ hiad respectively. 

The formula for the three-body decay is given by [1]: 


where E and p are the energies and momenta in the Iv 
rest frame. At the first equality we have used relations 
Pi = — p v and pv = Pd, while at the third equality 
we have used relations obtained by neglecting masses of 
leptons: 

E l =E v = \p v \ = -f!-. (16) 


r >= f m J dt 4 e i rD i 2 ' 

(ii) 

where P cm is the momentum of the final state vector 
meson in the D rest frame and p v is the momentum of the 
neutrino in the Iv rest frame, both of which are evaluated 
as 


P — 

1 cm — 


X 1/2 {m 2 Dl [M;^] 2 , m 2 v ) 


2m d 


( 12 ) 


Pv = 


A 1/2 ([M i ( ^ ) ] 2 , mf, ml) 


2 M 


(lls 


(13) 


with the Kallen function \(x, y, z) = x 2 + y 2 + z 2 — 
2 xy — 2 yz — 2 zx and the vector meson mass my. The 
tilde on characters for leptons indicates that they are 
evaluated in the Iv rest frame. The solid angles Q and 
f \ v are for the vector meson in the D rest frame and for 
the neutrino in the Iv rest frame, respectively, and 
is the Iv invariant mass. The integral range of is 

[nii+m„, mo — mv]- Substituting the squared amplitude 
with that in Eq. (10), we obtain 


r, = 


|g f | 2 


327 x b m 3 D mv 
2 


dM^P cm Pv dfl d£l v 


V, 


had 


{Pl • PD) {Pv ■PV )• 


(14) 


The energies and momentum of hadrons in the Iv rest 
frame can be exactly evaluated as 


E d = 


m2 D + [M^] 2 - my 


2 M; 


(Iv) 


(17) 


Ey = 


l D 


-[M! 


( il/ )l2 


2 M 


(Iv) 


(18) 


and \p D \ 2 = E'b — m 2 D . As a consequence, we have 
12 


r, = 


GfV<” 


87r 3 m| ) ?nv r 


dM^>P cm p v [M^] 2 


x Ie d E v - - \pd \ 2 


(19) 


where we have performed the integral of the solid angle 

n. 

In a similar way, we can evaluate the decay width for 
the four-body final state. The formula for the four-body 
decay is given by 



In general, the hadronization part Wad nray depend on where M^ t ' > is the invariant mass of the two-meson sys- 
the energy and scattering angles, and hence one cannot tern {hh), P'. m is the center-of-mass momentum of the 
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two-meson system in the D rest frame and ph is the mo¬ 
mentum of a meson in the hh rest frame, both of which 
are evaluated as 


P' 

cm 


A 1 / 2 (m 


2 

D ’ 


[aCi 2 . 

2m d 


( 21 ) 


Ph 


a i/2 m£ v h) r, m 


2 M: 


(hh) 



( 22 ) 


with the meson masses rrih and m' h . The momentum of 
the neutrino in the Iv rest frame p v is given in Eq. (13). 
The solid angles fl' and (lh are for the two-meson system 
in the D rest frame and for a meson in the hh rest frame, 
respectively. The tilde on characters for mesons indicates 
that they are evaluated in the hh rest frame. Since we 
are interested in the meson-meson invariant mass distri¬ 
butions for the semileptonic D decay, we calculate the 
differential decay width dTi/dM^^. Then in a similar 
manner to the case of the three-body decay, we have 


dT 4 




'jhh) / dM inv )P cmP^[ M i 


dM^ 32n 5 mpM- inv 


r(M 12 
J inv J 


x ( E d E s - -|pd | 2 ) , (23) 


where we have performed the integrals with respect to 
the solid angles Cl' and f lh- We mention that will 
be simply constructed as well, so that this can be put 
out of the integral, as we will see in the next subsection. 
The two-meson invariant mass can take a value 

within [ mh+m ' h , mo — mi — m v \, while the integral range 

of is [mi + m„, mo — The energies and 

momentum of hadrons in the parentheses can be exactly 
evaluated as 


structure. For the scalar mesons, we employ the chi¬ 
ral unitary approach [8-16], in which the scalar mesons 
are dynamically generated from the interaction of two 
pseudoscalar mesons governed by the chiral Lagrangians. 
Therefore, in this picture the light quark-antiquark pair 
after the W boson emission gets hadronized by adding 
an extra qq with the quantum number of the vacuum, 
uu + dd + ss, which results in two pseudoscalar mesons 
in the final state [see Fig. 1(b)]. Then the scalar mesons 
are obtained as a consequence of the final state interac¬ 
tion of the two pseudoscalar mesons as diagrammatically 
shown in Fig. 2. For the vector mesons, on the other 
hand, hadronization with an extra qq is unnecessary since 
they are expected to consist genuinely of a light quark 
antiquark pair [see Fig. 1(a)]. 


1. Scalar mesons 


First we consider processes with the scalar mesons in 
the final state as the dynamically generated resonances. 
The basic idea of the hadronization with an extra qq with 
the quantum number of the vacuum has already shown 
in Refs. [7, 17-19]. We start with the qq matrix M: 


M = 


uu ud us 
du dd ds 
y su sd ss 


(26) 


One can easily check that this matrix has the property 


M ■ M = M(uu + dd + ss). (27) 


With this property, a q/qp pair after the W boson emis¬ 
sion can be added by an extra qq to be 


dfdf -> {M ■ M) f p, (28) 


E d 


L D 


+ [ M «] 2 - 


2 M 


(lu) 


(24) 


where / denotes the flavor of light quarks: q\ = u, <72 = d, 
and <73 = s. Next we rewrite the matrix M in terms of 
the matrix <f> for pseudoscalar mesons 


Es 


m 2 D 


[M^] 2 - [M^l 

2 M<£> 


and \p D \ 2 = E 2 d - m 2 D 



B. Hadronizations 

Next we fix the mechanism for the appearance of the 
scalar and vector mesons in the final state of the semilep¬ 
tonic decay. We here note that, for the scalar and vec¬ 
tor mesons in the final state, the hadronization processes 
should be different from each other according to their 


FIG. 2: Diagrammatic representation of the direct plus rescat¬ 
tering processes for two pseudoscalar mesons. The solid and 
dashed lines denote quarks and pseudoscalar mesons, respec¬ 
tively. The shaded ellipses indicate the hadronization of a 
quark-antiquark pair into two pseudoscalar mesons, while 
the open circle indicates the rescattering of two pseudoscalar 
mesons. 
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75” 

7 r 

A'- 


i 

V6 77 


7T 1 


75-° 


-VS*!' 

K° 


7^ 


K+ 

K° 




¥) 


(29) 


where we have taken into account the rj-rf mixing in 
a standard way [44]. In this scheme we can calculate 
the weight of each pair of pseudoscalar mesons in the 
hadronization. Namely, the ss pair gets hadronized as 
ss(uu + dd + ss) = ((f) ■ (f>) 33, where 

(cf) ■ (f)) 33 = K~K+ + K°K° + (30) 

Here and in the following we omit the r\' contribution 
since v( is irrelevant to the description of the scalar 
mesons due to its large mass. In similar manners, the 
ds, sd , dd, su, and du pairs get hadronized as 

(</> ' ^) 23 = 7T _ A' + — -^=7T°A' 0 , (31) 

((j> ■ <j>) 32 =K~ 7 T+ - =K°n °, (32) 


1 


1 


(</> • </)22 = 7T 7T+ + ^7T U 7T U + ^T]T] - \l ^Tj + K°K° 


(33) 


and 


(/> • </)31 = ~^ 7T ° k + n K °' 


(< f >-< t >)21 = 77 + K°K , 


(34) 

(35) 


respectively. 

By using these weights, we can express the hadroniza¬ 
tion amplitude for the scalar mesons, in terms of 

two pseudoscalar mesons. For instance, we want to recon¬ 
struct /o(500) and / o (980) from the 7r + 7r _ system in the 
Df —> 7 r + 7 r l + vi decay. Because of the quark configu¬ 
ration in the parent particle D+, in this decay the 7 r + 7 r _ 
system should be obtained from the hadronization of the 
ss pair and the rescattering process for two pseudoscalar 
mesons, as seen in Fig. 2, with the weight in Eq. (30). 
Therefore, for the Df —> l + V( decay mode we can 

express the hadronization amplitude with a prefactor C 
and the Cabibbo-Kobayashi -Maskawa matrix elements 
V cs as 


^had^s"’ 7r+7r ] — CV cs ^G K + K -T K+K -_> n+n - 

+G k okoT k oko^.i I +-k- + g • 2 • -Gjjj• 

(36) 


In this equation, the decay mode is abbreviated as 
[D+, 7 T + 7 r ], and G and T are the loop function and 
scattering amplitude of two pseudoscalar mesons, respec¬ 
tively, whose formulation are given in Sec. IIC. We have 
introduced extra factors 2 and 1/2 for the identical par¬ 
ticles r/r/. The former factor 2 comes from the two ways 
of annihilating the 7777 operator in Eq. (30) by the I 7777 ) 
state as in the usual manner for effective Lagrangians, 
while the latter one 1/2 is the symmetry factor for the 
7777 loop. The scalar mesons /o(500) and /o(980) appear 
in the rescattering process and exist in the scattering am¬ 
plitude T for two pseudoscalar mesons. It is important 
that this is a Cabibbo favored process with V cs . Further¬ 
more, since the ss pair is hadronized, this is sensitive to 
the component of the strange quark in the scalar mesons. 
In this study we assume that C is a constant, and hence 
the hadronization amplitude V/7 is a function only of 
the invariant mass of two pseudoscalar mesons. Here we 
emphasize that the prefactor C should be common to all 
reactions for scalar meson productions, because in the 
hadronization the SU(3) flavor symmetry is reasonable, 
i.e., the light quark-antiquark pair q/cff hadronizes in 
the same way regardless of the quark flavor /. In this 
sense we obtain 


^hadl^s"’ K + K ] — CV cs ^1 + G k + k -T k + k -^ k + k - 


+G K o k oT k o K°^K+K + g ' 2 • ^ G r j r jT^_^j^+ 


(37) 


for the Df —> K + K l + Vi decay. In this case we have 
to take into account the direct production of the two 
pseudoscalar mesons without rescattering (the first dia¬ 
gram in Fig. 2), which results in the unity in the paren¬ 
theses. On the other hand, for the Df n~ K + l + vi 
decay mode the n~ K + system should be obtained from 
the hadronization of ds and hence this is a Cabibbo sup¬ 
pressed decay mode. The hadronization amplitude is ex¬ 
pressed as 


^hadl^s" ’ 71 K + ] — CVcd ^1 + G n -K+T- k- K+ -> 7 r- K+ 


--^Gt t O K oT 7i O K 0^t t - K + ) . 


(38) 


In a similar manner we can construct every hadroniza¬ 
tion amplitude for the scalar meson. The resulting ex- 



7 


pressions are as follows: 


V&[D + , t t + K~] = V^[D\ tt-K u ] = CV cs A„ k , (47) 


K 5 . 


V&[D+, 7T+7T-] = CV cd ( 1 + G n+n -T„ +v -^ + „- 


x • xG^o.-oT, 


1 1 


7T U 7T U -*■ 7T U 7T U —>-7T 


2 2 
+ G K °K°T K OK°->n+n- ) > 


7i T ^ • 2 • 2^ m '^ , ' nr, ^' K+lr ~ 

(39) 


Kal^ + , A] = CT c 


cd 


2 /2 

3 V 3^ 7 r ° , ^ ,r 0 l 7 _! ''"’ 0l 7 


^KORoTKOKO^nOri ! 


(40) 


^[£>+, A'+A-] = GV cd ^G 7 r+ 7 r -T OT+7r -^^ +K - 

+ 2 ' ^ ’ I^G^O^oT^o^o^r+r- + — • 2 • -G m T m ^ K + K - 


(41) 


^\l -^G 7 r o rl T 7 T o rl - > .K+K- + GroroT k o k o^r+r- J , 


^had[^ + ! W + K ] — CV ca ( 1 + G w + K~T-k+ K~ ->tt+ if- 


~~^=G n o K oT n o K o_ t . n +K- I i 


(42) 


L 7T T) —^7T T) 


V&[D°, n-r,} = CV cd (^= + -^=G n - v T„ 

k° K ~ r ^K°K~ r] ) , (43) 


v^ d [D°, at ° at ~] = cv cd (1 + ^G r ,r rHJf o K - 
+ G k o k -T k o k -^k«k- ) (44) 


Kad^°, 7T-X 0 ] = GV CS (1 + -4<W-T, 


V2 

Gtt-RoT^-K°->-k-K° 


t °K~ — Vtt K° 


(45) 


where A^k is a function of the invariant mass of two 
pseudoscalar mesons and is defined with the scattering 
amplitude in the isospin basis as 

AttK = 1 + G n kT^k(I=1/2)^hvK(I=1/2) • (48) 

In a similar manner, we simplify the hadroniza- 
tion amplitudes V^[.D + , n°rj\, Vhad[-D 0 , 7 r_? 7 ]i an( i 
V^[D°, K°K~] as 

Kad [D + , A] = --Lv^ d [D°, n- v ] = -)JlCV cd B nr ,, 

(49) 

Kad l D °i K ° K ~] = GV cd B KR , (50) 

with 

a/3 

-®7tt; = 1 + G^T^^^ ^“G ^ R Tr r (/=1)—(51) 

Brk =1 + GxK^/fK(/= 1 )_>ifK(j= 1 ) 

2 

~ ~^=G vri T nv _t K K( I -. 1 y (52) 

From the above expressions one can easily spec¬ 
ify Cabibbo favored and suppressed processes for the 
semileptonic decays into two pseudoscalar mesons, which 
are listed in Table II. 

Finally we note that the use of a constant G factor in 
our approach gets support from the work of Ref. [41]. The 
evaluation of the matrix elements in these processes is 
difficult and problematic. There are however some cases 
where the calculations can be kept under control. For 
the case of small recoil, namely when final pseudoscalars 
move slow, it can be explored in the heavy meson chiral 
perturbation theory [45]. Detailed calculations for the 
case of semileptonic decay are done in [41]. There one 
can see that for large values of the invariant mass of the 
lepton system the form factors can be calculated and the 
relevant ones in s wave that we need here are smooth in 
the range of the invariant masses of the pairs of mesons 
that we use here. To be able to use this behaviour we 
should prove that in our case the invariant masses of the 
lepton pair are large, but indeed, it was shown in the 
study of the semileptonic B decays [19] (and can be done 
also here) that the mass distribution of the lepton pair 
accumulates at the upper end of the phase space. There 
is also another limit, at large recoil, where an approach 
that combines both hard-scattering and low-energy in¬ 
teractions has been developed is also available [37], but 
this is not the case here. 


The hadronization amplitudes V] 1 ^[D+, 7r AT+], 
D h ^[D + , 7 t+A' _ ], and [D°, 7r _ A'°] are further 

simplified by using the isospin symmetry as 

V^[Dt, 7T-K+] = CV cd A„ K , (46) 


2. Vector mesons 

Next we consider processes with the vector mesons 
in the final state. As we have already mentioned, 



hadronization with an extra qq is unnecessary for the 
vector mesons. As a consequence, we can formulate the 
hadronization amplitude for vector mesons, in a 

very simple way. 

In order to see this, we consider the semileptonic de¬ 
cay Df 0(1020) l + vi as an example. The decay pro¬ 
cess is diagrammatically represented in Fig. 1(a), and the 
hadronization amplitude can be expressed with a 
prefactor C' and the Cabibbo-Kobayaslii -Maskawa ma¬ 
trix element V cs as 

Kad [Df, <t>\ = C'V cs , (53) 

where the decay mode is abbreviated as [D+, <j>] in the 
equation. Here we emphasize that the prefactor C' 
should be common to all reactions for vector meson pro¬ 
ductions, as in the case of the scalar meson productions, 
because the SU(3) flavor symmetry is reasonable in the 
hadronization, i.e., the light quark-antiquark pair q/qp 
hadronizes in the same way regardless of the quark flavor 
/. We further assume that C' is a constant again. This 
formulation is straightforwardly applied to other vector 
meson productions and we obtain the hadronization am¬ 
plitude for vector mesons: 


V^[Df, K*°] 

= C'V cd , 

(54) 

V^[D + , K*°\ 

| = CVcs, 

(55) 

v£\[d + , p° } = 

~T2 C ' Vcd ' 

(56) 

V^[D + , u] = 

7=2 C ' Vcd ’ 

(57) 

V&[D°, K*~] 

= -C'V cs , 

(58) 

Kal [D°, P~] 

= C'v cd , 

(59) 


where we have used K *, p and u> states in the isospin basis 
summarized in Appendix A. We note that these equa¬ 
tions clearly indicate Cabibbo favored and suppressed 
processes with the Cabibbo-Kobayaslri-Maskawa matrix 
elements V cs and V cd , respectively. 


where i, j, and k are channel indices, s is the Mandelstam 
variable of the scattering, V is the interaction kernel, and 
G is the two-body loop function. For the hadronization in 
the previous subsection we need three types of coupled- 
channels systems: the ( Q , S) = (0, 0) system, for which 
we introduce six channels labeled by the indices i = 1 , 
..., 6 in the order 7r + 7r _ , 7r°7r°, K + K~, K°K°, ryri, and 
7 r° 77 , the KK(I = l)-7 rrj system, and the nK(I = 1/2)- 
r/K system. 

In this study the interaction kernel V.- t j = Vj, is taken 
as the simplest one, that is, the leading-order s-wave in¬ 
teraction obtained from the chiral perturbation theory. 
The interaction kernel for (Q, S) = (0, 0) is summarized 
as 


Vn = 2 Vi 3 = 2Hi4 = 2V2V 23 = 2s/2V 2A 

= ^33 = 21^34 = V44 = — 7^2’ 


Vi 2 

Vl 5 

V 16 

V 35 


S ~~ 

~ V2p ’ 

_ 1 

-jj-V22 = V2V25 = y=V66 = 

V 26 = V 56 = 0 , 

9s — 2 ml — 6 rrii 

--SW 1 



V 36 

V 55 


~v 46 = - 


9s — m%. — 8 m 2 K — 3 m 2 


12 v / 3/ 2 


7 m 2 — 16 m 2 K 

W 2 : 


(61) 


where / is the pion decay constant. One must remem¬ 
ber that in the chiral unitary approach when calculat¬ 
ing T = (1 — VG)~ X V one uses the unitary normaliza¬ 
tion (l/v / 2 )| 7 r° 7 r°) and (l /\f2)\r)rj) for identical particles, 
which allows to use the general formula in coupled chan¬ 
nels. At the end the good normalization of the external 
particles must be restored and these are the amplitudes 
that appear to Eq. (38) and following ones. 

For the KK(I = scattering, the interaction ker¬ 

nel can be written in terms of the interaction kernel for 
the (Q, S ) = (0, 0) system shown above (see the isospin 
basis summarized in Appendix A): 


Vkk(i=i)->kk(i=i) — 2 (^33 — 2 F 34 + V44) , (62) 


C. Scattering amplitudes of two pseudoscalar 
mesons in chiral unitary approach 


^ KK(I=l)-n-nri — ~ (^36 — ^46) j 


(63) 


For the scattering amplitude of two pseudoscalar 
mesons, we employ the so-called chiral unitary ap¬ 
proach [8-16], which we briefly explain in this subsection. 
In this approach we solve a coupled-channels Betlre- 
Salpeter equation in an algebraic form 

Tij(a) = Vij(s) + Y, V ik (s)G k (s)T kj (s), (60) 

k 


V-TT^—HTr] — 1*66 • (64) 

For the ttK(I = 1/2 )-pK scattering, the interaction ker¬ 
nel is expressed as 

VnK(I=l/2)^nK(I=l/2) = g^p [ — ® s2 + ^( m n + TO I:) S 

+ 3 (m n - m K ) 2 ] , 


(65) 
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V-kK(I=1/2)<+tiK 


~ 04s f 2 ^ s2 2m tf + 3m 2 )s 

—9 {ml- m 2 K )(m 2 K - m 2 )\ , ( 66 ) 


Vr,K^r,K t 9s2 + 2 ( 2m2 “ 9m K ~ 3m ^) S 

+ 9 (m K ~ m.r,) 2 ] . ( 67 ) 

For the loop function G, on the other hand, we use the 
following expression: 

r ( . = . [ d 4 q _ 1 _ 1 _ 

* 1 J (27r) 4 q 2 — m 2 + *0 (P — q ) 2 — mf + i 0 

( 68 ) 

where = (y/s, 0) and m.; and m\ are masses of pseu¬ 
doscalar mesons in channel *. In this study we employ a 
three-dimensional cut-off g max as 


r ( x f d 3 q U}j(q) + ui'ijq) _ 6 (q max \q\) 

J (27r) 3 2 Ui(q)cj-(q) s - [usi(q) + v'(q )] 2 + i O’ 

(69) 

In this expression we have performed the q° integral and 
coi(q) = yjrrif + q 2 and u)'(q) = yjm f + q 2 are the on- 
shell energies. 

In this framework, with a small number of free pa¬ 
rameters we can reproduce experimental observables of 
meson-meson scatterings fairly well. In this study we 
take the model parameters of the chiral unitary approach 
as / = 93 MeV and q max = 600 MeV, which dynamically 
generates resonance poles in the complex energy plane: 
453 - 253* MeV for / o (500), 982 - 5* MeV for / o (980), 
and 721 — 236* MeV for /Vg(800). The ao(980) appears 
as a cusp at the KK threshold. 



FIG. 3: Meson-meson invariant mass distributions for the 
semileptonic decay Df —> PPe + v e with PP = n + n~, 
K + K~, and K + in s wave. We multiply the ty~ K + mass 
distribution, which is a Cabibbo suppressed process, by 10. 



III. NUMERICAL RESULTS 

Now let us calculate the semileptonic decay widths of 
D mesons into scalar and vector mesons. As we have 
formulated, we have only one model parameter for scalar 
and vector meson productions, respectively. Namely one 
can calculate the decay widths of the scalar meson pro¬ 
ductions with one common parameter (7, and similarly 
C’ for the vector meson productions. 

First we consider the scalar meson production in 
Sec. Ill A, and then move to the vector meson production 
in Sec. IIIB. Finally in Sec. IIIC we compare the two 
contributions of the mass distributions from the scalar 
and vector mesons. 


A. Production of scalar mesons 

In order to calculate the branching fractions of the 
scalar meson productions, we first fix the prefactor con¬ 
stant C so as to reproduce the experimental branching 
fraction which has the smallest experimental error for the 


FIG. 4: Meson-meson invariant mass distributions for the 
semileptonic decay D + —> PPe + v e with PP = 7r + 7r _ , n°r /, 
K + K~, and tv + K~ in s wave. We multiply the 7r + 7r _ , 7 t °r), 
and K+K~ mass distributions, which are Cabibbo suppressed 
processes, by 10. 


process with the s-wave two pseudoscalar mesons, that is, 
B[D+ -> (t r+Ar-) s . wave e+i/ e ] = (2.32 ±0.10) x 10" 3 . By 
integrating the differential decay width, or mass distri¬ 
bution, dT 4 /dM ^' 1 in an appropriate range, in the case 
of 7 t + A'~ [rn.Tr + rriK, 1 GeV], we find that C = 4.597 can 
reproduce the branching fraction of (tt + K~) s . wave e + u e . 

By using the common prefactor C = 4.597, we can cal¬ 
culate the mass distributions of two pseudoscalar mesons 
in s wave for all scalar meson modes, which are plotted 
in Figs. 3, 4, and 5 for Df, D + , and D° semileptonic de¬ 
cays, respectively. We show the mass distributions with 
the lepton flavor l = e; the contribution from l = /j is 
almost the same as that from l = e in each meson-meson 
mode due to the small lepton masses. In each figure we 
multiply the mass distributions which are Cabibbo sup- 
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FIG. 5: Meson-meson invariant mass distributions for the 
semileptonic decay D° —> PPe + v e with PP = ir~ri, K°K ~, 
and 7 t~K° in s wave. We multiply the -K~r) and K°K~ mass 
distributions, which are Cabibbo suppressed processes, by 10. 


pressed processes by 10 so that we can easily compare 
the shape of the mass distributions. As one can see, the 
largest value of the mass distribution dTi/dM^^ is ob¬ 
tained in the Df —> TT + ir~e + i , e process, in which we can 
see a clear / o (980) peak. It is interesting to note that 
in the Df n + TT~e + v e process we find a clear / 0 (980) 
signal while the /o(500) contribution is negligible, which 
strongly indicates a substantial fraction of the strange 
quarks in the /o(980) meson, as we will discuss later. For 
the Df semileptonic decay we also observe a rapid en¬ 
hancement of the I\ + K~ mass distribution at the thresh¬ 
old, as a tail of the /o(980) contribution, although its 
height is much smaller than the 7 r + 7 W peak. For the D + 
and D° semileptonic decays, we can see the tt + K~ and 
7T ~K° as Cabibbo favored processes, respectively. We 
note that the n + K~ and n~K° mass distributions are 
almost the same due to isospin symmetry. It is interest¬ 
ing to see that the shape of the 7r + A'“ and n~ K° mass 
distributions is determined by, in addition to the A"q( 800) 
resonance, the kinetic factor of the squared decay ampli¬ 
tude. Namely, we have the matrix element of Eq. (10) 
that is roughly proportional to | p u \ 2 and this momentum 
gets bigger the smaller the meson-meson invariant mass. 
This kinetic factor of the squared decay amplitude affects 
the 7 r + 7 r - distribution in the D + semileptonic decay in a 
similar manner, and also provides more weight at low in¬ 
variant masses for the shape for mj in Figs. 4 and 5 than 
the n°r] distributions in the D° K 0 tt°t] decay evalu¬ 
ated in Ref. [18]. The nr] mass distributions in Figs. 4 and 
5 of the D + and D° decays show peaks corresponding to 
ao(980), but its peak is not high compared to the /o(980) 
peak in the 7r + 7r~ mass distribution of the Df decay since 
they are obtained in Cabibbo suppressed processes. The 
D + —» 7T + 7r~e + i' e decay is Cabibbo suppressed and it 
has a large contribution from the /o(500) formation and 
a small one of the /o(980), similar to what is found in 



FIG. 6: 7r + 7r _ invariant mass distribution for the semilep¬ 
tonic decay D+ —> 7r + 7r _ e + i' e . The theoretical calculation is 
folded with the size of experimental bins, 25 MeV. The exper¬ 
imental data are taken from Ref. [27] and are scaled so that 
the fitted Breit-Wigner distribution (dashed line) reproduces 
the branching fraction of B[D+ — >• /o(980)e + Z7 e , / o (980) —► 
7T + 7r _ ] = 0.2% by the Particle Data Group (see Table I). 

the B° —> J/i/)ir + Tr~ decay in Ref. [7]. A different way to 
account for the PP distribution is by means of dispersion 
relations, as used in Ref. [41] in the semileptonic decay 
of R, where the 7T + 7r” s-wave distribution has a shape 
similar to ours. 

The theoretical 7r + 7r _ mass distribution of the semilep¬ 
tonic decay D s —> e + v e is compared with the 

experimental data [27] in Fig. 6. We note that we 
plot the figure in unit of ns _1 /GeV, not in arbitrary 
units. The theoretical mass distribution is folded with 
25 MeV spans since the experimental data are collected 
in bins of 25 MeV. The experimental data, on the other 
hand, are scaled so that the fitted Breit-Wigner distri¬ 
bution reproduces the branching fraction of B[Df —> 
/o(980)e + ^ e , / o (980) —> 7r+7r~] = 0.2% [1], The mass 
and width of the Breit-Wigner distribution are fixed as 
M = 966 MeV and F = 89 MeV, respectively, taken 
from Ref. [27]. In Fig. 6 we can see a qualitative cor¬ 
respondence between the theoretical and experimental 
signals of /o(980). We emphasize that, both in exper¬ 
imental and theoretical results, the 7r + 7r - mass distri¬ 
bution shows a clear / 0 ( 980) signal while the /o(500) 
contribution is negligible. This strongly indicates that 
the /o(980) has a substantial fraction of the strange 
quarks while the /o(500) has a negligible strange quark 
component. It is interesting to recall that the appear¬ 
ance of the /o(980) in the case one has a hadronized 
ss component at the end, and no signal of the /o(500), 
is also observed in R° and B° decays in Refs. [2-6]. 
The explanation of this feature along the lines used in 
the present work was given in Ref. [7]. However, al¬ 
though the peak height of the /o(980) is very similar, the 
Breit-Wigner fit would provide larger branching fraction 
B[Df —> /o(980)e + r' e ] than the theoretical one. Actually, 
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integrating the theoretical mass distribution in the range 
[0.9 GeV, 1.0 GeV], we obtain the branching fraction 
B[D+ -> / o (980)e+i/ e ; / o (980) tt+tt"] = 5.10 x 10“ 4 , 
which is about four times smaller than the experimental 
value 2.00 x 10~ 3 . Actually, in the experimental analysis 
of Ref. [27] different sources of background are considered 
that make up for the lower mass region of the distribu¬ 
tion. The width of the /o(980) extracted in the analysis 
of Ref. [27] is T = 91 433 A 3 MeV, which is large com¬ 
pared to most experiments [1], including the LHCb ex¬ 
periment of [46], although the admitted uncertainties are 
also large. One should also take into account that, while 
a Breit-Wigner distribution for the / 0 (980) is used in the 
analysis of Ref. [27], the large coupling of the resonance 
to KK requires a Flatte form that brings down fast the 
7 r + 7 r~ mass distribution above the KK threshold. Our 
normalization in Fig. 6 is done with the central value of 
the B[D + —> (tt + K - ) s _ wave e + z/ e ] and no extra uncertain¬ 
ties from this branching fraction are considered. Yet, we 
find instructive to do an exercise, adding to our results 
a “background” of 10 ns -1 /GeV from different sources 
that our calculation does not take into account, and then 
our signal for the /o(980) has a good agreement with the 
peak of the experimental distribution. 

As mentioned above, the value extracted in [27] for 
the / 0 (980) signal is tied to the assumptions made, in¬ 
cluding parts of the background that lead to a very 
large width of the resonance, assuming a Breit-Wigner 
shape, etc. Actually, in a more recent paper [47] the 
same CLEO data of [27] are reanalyzed taking a band of 
/o(980) masses within 60 MeV of 980 MeV and assuming 
a Flatte form of the resonance and a rate for £>[D+ —> 
/o(980)e+z4, / o (980) -> tttt] = (0.13 ± 0.02 ± 0.01)% is 
obtained. This value is about a factor of two smaller than 
the one reported in [27] and more in agreement with our 
results. 

Next we consider the differential decay width with re¬ 
spect to the squared momentum transfer q 2 , which coin¬ 
cides with the squared invariant mass of the lepton pair: 
q 2 = [M^] 2 . The differential decay width for the scalar 
meson production is expressed as 


dr 4 


|g f | 2 


647T 5 7 


dM 


( hh ) 


K 


M 


had 


P'phPuM 


l D 


Ml 

x (e d E s - -|pd | 2 


(hh) 


(70) 


This differential decay width was experimentally ob¬ 
served in Ref. [27] for the Df / o (980)e + iz e decay mode 
followed by / 0 (980) —> ir + ir~. In this study we compare 
our theoretical value for this decay mode with the ex¬ 
perimental data in Fig. 7. The range of the integral for 
is [0.9 GeV, 1.0 GeV]. As one can see, we can to 
some extent reproduce the shape of the differential decay 
width dT^/dq 2 in experiment, but the absolute value of 
the theoretical calculation is several times smaller than 
the experimental one. This can be, as we have explained, 
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FIG. 7: Differential decay width of the D/ —r /o(980)e + i' e 
decay mode followed by /o(980) — »• 7r + 7r _ , with q 2 = [M^] 2 . 
The experimental data are taken from Ref. [27]. The exper¬ 
imental points should be rescaled dividing by about a factor 
of two if the absolute rate for the /o(980) production of the 
reanalysis of Ref. [47] were used. 


solved by introducing background contributions when ex¬ 
tracting the amount of the ,/b(980) signal from experi¬ 
mental data. Actually, as we have commented before, 
the reanalysis of [47] leads to absolute values of the rate 
for the /o(980) production about a factor of two smaller, 
and again if we scale the q 2 distribution of in Fig. 7 by 
the factor the agreement is much better. 

Moreover, integrating the mass distributions we cal- 


TABLE III: Branching fractions of semileptonic D decays into 
two pseudoscalar mesons in s wave. The branching fraction 
of the D + —> (7r + K~ )s- W a Ve e + v £ mode is used as an input. 


Mode 

Df 

Range of [GeV] 

l = e 

l = /j, 

7T + 7r _ 

K+K~ 

tv~I< + 

[0.9, 1.0] 

[2mx, 1.2] 

[m,r + rtiK, 1.0] 

5.10 x 10^ 4 
1.42 x 10^ 4 

8.11 x 10^ 5 

4.71 x 10^ 4 
1.30 x 10^ 4 
7.63 x 10^ 5 

Mode 

D+ 

Range of M^ [GeV] 

l = e 

l = n 

7T + 7r _ 

7T°?7 

K+K~ 

7 r + K~ 

[2rrj.Tr, 1.0] 

[m,, + m n , 1.1] 
[2m*:, 1.2] 

[m,r + m,K, 1.0] 

5.11 x 10^ 4 
6.37 x 10^ 5 
2.24 x 10^ 6 
2.32 x 10^ 3 

4.85 x 10~ 4 

5.86 x 10^ 5 
2.01 x 10^ 6 
2.16 x 10~ 3 

Mode 

D° 

Range of M^ [GeV] 

l = e 

l = n 

7T _ 77 

K°K~ 

7 t~K° 

[m-n + m n , 1 . 1 ] 
[2m*:, 1.2] 

[m,, + m K , 1.0] 

4.93 x 10” 5 
5.47 x 10^ 6 
8.99 x 10^ 4 

4.53 x 10~ 5 
4.88 x 10^ 6 
8.38 x 10^ 4 
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TABLE IV: Branching fractions of semileptonic D decays into 
vector mesons. 


Mode 

D+ 

l = e 

l — p 

0 (1020) 

2.12 x 10 -2 

1.94 x 10“ 2 

K* (892)° 

2.02 x 10” 3 

1.89 x 10~ 3 


D+ 


Mode 

l = e 

l — fi 

K*( 892)° 

5.56 x 10~ 2 

5.12 x 10” 2 

p(770)° 

2.54 x 10~ 3 

2.37 x 10” 3 

w(782) 

2.46 x 10~ 3 

2.29 x 10~ 3 


D° 


Mode 

l = e 

l = p 

K* (892)' 

2.15 x 10~ 2 

1.98 x 10” 2 

P(770)- 

1.97 x 10~ 3 

1.84 x 10~ 3 


culate the branching fractions of the semileptonic D 
mesons into two pseudoscalar mesons in s wave, which 
are listed in Table III. We note that the branching frac¬ 
tion B[D + — > (n + K _ ) s . wave e + ^ e ] = 2.32 x 10 ~ 3 is used 
as an input to fix the common constant, C = 4.597. 
Among the listed values, we can compare the theoret¬ 
ical and experimental values of the branching fraction 
B[D+ (K + K~) s . wave e + p e ]. Namely, in Ref. [24] this 
branching fraction is obtained as (0.22+o;Jg ± 0.03)% of 
the total Z)+ —> K + K~e + v e , which is dominated by the 
0(1020) vector meson. This indicates, together with the 
branching fraction Df —> 0 ( 1020 )e + is e , we can estimate 
B[D+ ->• {K+K~) sm e+v e ] = (5.5t 3 ;}) x 10“ 5 . Theo¬ 
retically this is 1.42 x 10 -4 . Although our value overesti¬ 
mates the mean value of the experimental data, it is still 
in 3cr errors of the experimental value. 


B. Production of vector mesons 

Let us move to the vector meson productions in the 
semileptonic D decays. For the vector mesons we fix the 
common prefactor C' so as to reproduce the 10 available 
experimental branching fractions listed in Table I. From 
the best fit we obtain the value C = 1.563 GeV, which 
gives X 2 /-/V<j.o.f. = 22.8/9 « 2.53. The theoretical values 
of the branching fractions are listed in Table IV and are 
compared with the experimental data in Fig. 8 , where we 
plot the ratio of the experimental to theoretical branch¬ 
ing fractions. We calculate the experimental branching 
fraction of the D + K(892)°l + vi (l = e and fi) process 
by dividing the value in Table I by the branching fraction 
£>[A"*(892) 0 —> K~n + ] = 2/3, which is obtained with 
isospin symmetry. As one can see from Fig. 8 , the exper¬ 
imental values are reproduced well solely by the model 
parameter C' with X 2 /-^d.o.f. ~ 2.53. 
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FIG. 8: Ratio of the experimental to theoretical branching 
fractions for the semileptonic D decays into vector mesons. 


Next for the —> 0(lO2O)e + z/ e decay mode we con¬ 

sider the differential decay width with respect to the 
squared momentum transfer q 2 , which coincides with the 
squared invariant mass of the lepton pair: q 2 = [M^] 2 . 
This differential decay width was already measured in an 
experiment [27] for the Df 0(1020 )e + v e decay mode. 
In a similar manner to the previous case, the differential 
decay width for the vector meson production is expressed 
as 


d£s 
dq 2 


C' t/M 

^h a d (jj.\ 

___ p f, My * 

1 §ir 3 m z D mv cm v lnv 


^ Ed By 



(71) 

In Fig. 9 we compare our result for this reaction with 
the experimental one. As one can see, our theoretical 
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FIG. 9: Differential decay width of the Df —>• 0(lO2O)e + i' e 
decay mode followed by 0(1020) — > K + K~, with q 2 = 
[A//// 2 . The experimental data are taken from Ref. [27]. The 
theoretical value is multiplied by the branching fraction of the 
0(1020) meson to K+K~, S[0(1O2O) -»• K+K~] = 48.9% [1]. 
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result reproduces the experimental value of the differen¬ 
tial decay width quantitatively well. This means that our 
method to calculate the semileptonic decays of D mesons 
is good enough to describe the decays into vector mesons. 

In this study we have not evaluated the D+ —> 
<j>( 1020)e + i/ e decay. This decay proceeds like the D+ —> 
w(782)e + i/ e decay that we have evaluated and one has a 
dd at the end. Since the <j> is ss then this is forbidden 
in our approach, at the tree level that we have consid¬ 
ered for the vector production. Experimentally, this rate 
is < 9 x 10 -5 . This is an upper bound about 30 times 
smaller than the rate of the omega production that we 
have evaluated. We do not want to go beyond, but can 
give some idea on how a finite rate could be obtained in 
our approach. For this one would have to hadronize the 
dd into a K°K°, then have a loop for K°K° propagation 
in p -wave and finally have the K°K° couple to the (f). 
Some technical details could be borrowed from the study 
of 4> —> 7T7T decay studied in [48] but one can get an indica¬ 
tion that the rate should be rather small by simply noting 
that the hadronization to meson-meson pairs has a re¬ 
duction factor, as one can see by comparing for instance 
/o(500) production with p production [49]. On the other 
hand, the coupling of <f> to K°K° is intrinsically small, as 
one can see from the 1.5 MeV partial decay width of this 
channel [comparatively the A(1232) partial decay width 
to the 7 tN channel would be about 15 MeV for a pion 
with the same momentum as the kaon in the (f) decay]. 
There are other factors to consider, but this can give us 
a feeling that the rate could be some orders of magnitude 
smaller than for omega production. 


C. Comparison between scalar and vector meson 
contributions 

Finally we compare the mass distributions of the two 
pseudoscalar mesons in s- and p -wave contributions. In 
the present approach the s-wave part comes from the 
rescattering of two pseudoscalar mesons including the 
scalar meson contribution, while the p -wave one from the 
decay of a vector meson. In this study we consider three 
decay modes: D+ —> n + ir~e + v e , D+ —> I\+ K~ e+v e , 
and D+ tx+K~ e+v e . The D+ —» 7r + 7r - e + i' e de¬ 

cay mode would have a large p-wave contribution from 
p( 770), but we do not consider this decay mode since it 
is a Cabibbo suppressed process. 

First we consider the D+ —> tt+t{- e+v e decay mode. 
This is a specially clean mode, since it does not have 
vector meson contributions and is dominated by the s- 
wave part. Namely, while the 7r + 7r - can come from a 
scalar meson, the primary quark-antiquark pair in the 
semileptonic D+ decay is ss, which is isospin 1 = 0 and 
hence the p(770) cannot contribute to the 7r + 7r~ mass 
distribution. The primary ss can be 0(1020), but it de¬ 
cays dominantly to KK and the 0(1020) — > n + ir~ decay 
is negligible. This fact enables us to observe the scalar 
meson peak without a contamination from vector meson 



3500 

j-, 

> 

0 

3000 


2500 



2000 


1500 

§ 

1000 

_ -"t 


500 

0 


0.95 



s wave 
p wave 


1.05 

M^ h> [GeV] 


1.1 


1.15 1.2 


FIG. 10: K+K invariant mass distribution for the semilep¬ 
tonic decay D+ —> K+K~ e+v e both in s and p waves. 


decays and discuss the quark configuration in the /o(980) 
resonance as in Sec. Ill A. 

Next let us consider the D+ K+K~e+v e decay 
mode. As we have seen, the K + K~ mass distribution 
in s wave is a consequence of the / 0 (980) tail. However, 
its contribution should be largely contaminated by the 
0(1020) -> K+K~ in p wave, which has a larger branch¬ 
ing fraction than the (K+K~) s . wave in the semileptonic 
decay. In order to see this, we calculate the p -wave 
K+K~ mass distribution for D+ —> I\+ K~ e+v e , which 
can be formulated as 


dT 3 2m Vt T 3 xfflF->M] 

-=-ini - 

dM^ " [M^Y-ml + imvVviM^y 

(72) 

where my is the vector meson mass and the energy de¬ 
pendent decay width Ty(M^^) is defined as 


r v {M^) = T v 


(73) 


P° S (M^ h) ) = 


A 1/2 ([M i ( “ ) ] 2 , ml, m'fi) 


2 Ml 


( hh ) 


(74) 


p on = AV2 K, m h) (75) 

2 my 

For the 0(1020) meson we take (0 = 4.27 MeV and 
B[</> —> K+K~] = 0.489 [1]. The numerical result 
for the (K+K~) p _ wave mass distribution is shown in 
Fig. 10 together with the {K+ A" _ )s-wave- From the 
figure, compared to the (K+ A' _ ) p - waV e contribution we 
cannot find any significant (K+K~) s _ wave contribution, 
which was already noted in the experimental mass dis¬ 
tribution in Ref. [24]. Nevertheless, we emphasize that 
the (K+K~) s . wave fraction of the semileptonic D+ de¬ 
cay is large enough to be extracted [24], Actually in 
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FIG. 11: 7 t + A'~ invariant mass distribution for the sernilep- 
tonic decay D + —> n + K~ e + v e both in s and p waves. 


Ref. [24] they extracted the {K + K - ) s - wave fraction by 
analysing the interference between the s- and p -wave 
contributions. This fact, and the qualitative reproduc¬ 
tion of the branching fractions in our model, implies that 
the /o(980) resonance couples to the KK channel with 
a certain strength, which can be translated into the I\K 
component in /o(980), in a similar manner to the KD 
component in £>* 0 (2317) [19, 50], in terms of the compos¬ 
iteness [51]. Anyway, in order to conclude the structure 
of the /o(980) more clearly, it is important to reduce the 
experimental errors on the (K + K~) s . W!LVe . 

Finally we consider the D + —> n + K~e + v e decay mode. 
In this mode the ( 7 r + A" _ ) s _ wave from the A"q( 800) and the 
( 7 r + A' - ) p _ wave from the A"*(892) are competing with each 
other. In a similar manner to the D+ —> K + K~e + v e 
case, we calculate the mass distribution also for the 
p -wave 7 T+AT - contribution dT 3 /dM with f^. = 
49.1 MeV [1], and the result is shown in Fig. 11. As 
one can see, thanks to the width of ~ 50 MeV for the 
AT* (892), the s-wave component can dominate the mass 
distribution below 0.8 GeV. We note that we would ob¬ 
tain an almost similar result for the D° n~K°e + u e 
decay mode due to isospin symmetry. 

As to the theoretical uncertainties, we can play a bit 
with the cut-offs used to regularize the loops, such that 
the masses of the states do not change appreciably. This 
exercise has been done a number of times and given us 
the feeling that within our models the uncertainties are 
below 10%. For the case of scalar production where we 
have a range of invariant masses and rely upon a constant 
production vertex C, the changes with the invariant mass 
in the primary form factors, prior to the final state inter¬ 
action of the mesons, as found in [41], would add some 
extra uncertainty. In total it would be fair to accept 
about 20 % uncertainties in this case in the limited range 
of energies that we move. 


IV. CONCLUSION 

In this study we have discussed the semileptonic de¬ 
cays of D mesons into light scalar and vector mesons. 
For the scalar meson production, we have formulated the 
semileptonic decay as the combination of two parts. One 
is the weak decay of the charm quark and the emission 
of a lepton pair via the W boson. The other is a sim¬ 
ple hadronization of light qq pair plus an extra qq from 
vacuum into two pseudoscalar mesons after the W boson 
emission, so as to generate the scalar mesons dynami¬ 
cally in the meson-meson final state interaction. The 
hadronization naturally gives the weight of each pair of 
pseudoscalar mesons in the decay process, which governs 
which scalar meson appears in the decay mode. For the 
vector mesons, on the other hand, we have not consid¬ 
ered the hadronization with an extra qq and have directly 
used the light qq pair after the W boson emission as a 
weight for the vector mesons, which are expected to be 
genuinely qq states. We note that we can specify flavors 
of quarks contained in the final state scalar and vector 
mesons by considering Cabibbo favored and suppressed 
processes. In addition, since the leptons interact only 
weakly, the semileptonic decay of the heavy meson to 
two light mesons +l + vi brings a suitable condition to 
measure effects of the final state interaction of the two 
light mesons. 

In our model of the semileptonic decay, the produc¬ 
tion yields of the scalar and vector mesons are respec¬ 
tively determined solely by constant prefactors C and C' 
as model parameters. Fixing C from the branching frac¬ 
tion of the D + —> ( 7 r + A" - ) s _ wave e + r 4 decay, we have cal¬ 
culated branching fractions of scalar meson productions. 
We have qualitatively reproduced the experimental value 
of the branching fractions of Df —» ( 7 T + 7 r _ ) s - wa v e e + t' e 
and (A' + A' _ ) s _ wave e + i' e decay modes. Some deviations of 
these branching fractions compared to the experimental 
values can be explained by taking into account the back¬ 
ground of the mass distribution for the 7 T + 7 r _ case and by 
the large experimental error for the K + K~ case. For the 
vector mesons, we have determined the constant C' so as 
to fit our numerical values to the available experimental 
values of the branching fractions, and we have reproduced 
the experimental values at a quantitative level. We also 
compared the mass distributions of the two pseudoscalar 
mesons in s- and p -wave contributions, which come from 
decays of the scalar and vector mesons, respectively. 

We have found that the Cabibbo favored decay mode 
Df f o {980)l+vi followed by /o(980) —» 7r + 7r~ and 
K + K~ is of special interest. For the / o (980) —> 
mode, we have found that there is no p -wave contami¬ 
nation from p(770) decay and hence it should be dom¬ 
inated by the s-wave part. Then, we have confirmed 
the experimental fact that the 7 r + 7 U mass distribution 
shows a clear / o (980) signal while the /o(500) contribu¬ 
tion is negligible. This strongly indicates that the /o(980) 
has a substantial fraction of the strange quarks while the 
/o(500) has a negligible strange quark component. For 
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the /o(980) —> K + K~ mode, on the other hand, the 
{K + K _ ) s . wave contribution is highly contaminated by 
the 0(1020) —> K + K~ decay in p wave. Nevertheless, 
the (K + K~) s . wave fraction of the semileptonic D+ decay 
is large enough to be extracted experimentally, which im¬ 
plies that the /o(980) resonance couples to the KK chan¬ 
nel with a certain strength and hence implies a certain 
amount of the KK component in /o(980). 
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with u = Y 7 ° and v = t^y 0 , and hence we have 
^u(p, s)u(p, s) = 

' (A2) 

^2v(p, s)v{p , s) = 

S 

where m is the mass of the field, j = 7 with y p being 
the Dirac 7 matrices, and + m 2 , p'j is the 

on-shell four-momentum of the solution. 

The identities for the Dirac matrices used in this study 
are summarized as follows: 

7°(7 M ) 1 7° = Y, (75) t = 75, (A3) 

tr [YYYY] = 4 (<T + g^g^), (A4) 

tr [y 5 YYYY] = -4 ie^ pa , (A5) 

tr [YYY] = tr [jnYYY] = 0, (A 6 ) 

where 75 = ij°YYY and e A ‘ 1/pcr is the Levi-Civita sym¬ 
bol with the normalization e 0123 = 1 . The Levi-Civita 
symbol satisfies the following identity 

e a ^e a0pa = -2 (<fX ~ «)• (A7) 

3. Isospin basis 


Appendix A: Conventions 

In this Appendix we summarize conventions used in 
this study. 


1. Metric and Lorentz indices 

In this article the metric in four-dimensional 
Minkowski space is g^ = g^ v = diag(l, — 1, — 1, — 1) 
and the Einstein summation convention is used unless 
explicitly mentioned. The scalar product of two vectors 
a p and 6 P is represented as a ■ b = = a°b° — a ■ b. 


2. Dirac spinors and matrices 

As the positive and negative energy solutions of the 
Dirac equation, we express the Dirac spinors respectively 
as u(p , s) and v(p, s), where p is three-momentum of 
the field and s represents its spin. The Dirac spinors are 
normalized as follows: 

v(p, s)v(p, s') = - 6 SS >, (Al) 


In terms of the isospin states 1 1, I 3 ), the phase conven¬ 
tion for pseudoscalar mesons is given by 

k + > = —11, 1), |A - } = 11/2, —1/2), (A 8 ) 

while other pseudoscalar mesons are represented without 
phase factors. As a result, we can translate the physi¬ 
cal two-pseudoscalar meson states into the isospin basis, 
which we specify as (J, J 3 ), as 


\KK(0, 0)) = --L\K+K-) - -J=| K°K°), 

(A9) 

|r?r?( 0 , 0 )) = | rm), 

(A10) 

\KK(1, 0)) = -±=\K+K~) + ±=\K°K°), 

(All) 

\KK(l, -1)) = -| K°K~), 

(A12) 

M(l, 0)) = | 7 r°? 7 >. 

(A13) 

1 ^( 4 , “!)) = \n~v), 

(A14) 


u(p, s)u(p , s') = 6 as >, 
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W 1/2, -1/2)) = -^=\tt°K°) - \n-K+), (A15) 

\nK(l/2, 1 / 2 )) = \[\U + K~) - (A16) 

|ttA'(1/2, — 1 / 2 )) = --L\n°K~) - \n~K°), (A17) 

Furthermore, the vector meson states are represented 
in terms of quarks as 


where V cq is the Cabibbo-Kobayashi-Maskawa matrix 
elements for the transition from the charm to light quark 
q. The W boson propagator with four-momentum is 
written as 

iP ^ ip) = p^-M, 2, + iO- (A23) 

with the mass of the W boson M w . The coupling con¬ 
stant gw and the mass of the W boson Mw are related 
to the Fermi coupling constant Gp as 


I P°) 


s/T 



I P ) = I du), 


(A18) 


G f 


9w 

4 y/2M*y' 


(A24) 


|w) 


-^\ uu ) + ^ dd ^ 


(A19) 


5. Physical constants 


|A"* 0 ) = |ds), |AT* 0 ) = |sd), | K*~) = -|su). (A20) 

4. Feynman rules 

The Wul coupling is expressed as 

~ iV ^i =i % (A21) 

with gw being the coupling constant of the weak inter¬ 
action, and the Wcq coupling as 

-iV£„ = ^ 7* V/L (A 22 ) 


In this article we use the following values for phys¬ 
ical constants. The Fermi coupling constant: Gp « 
1.166 x 10 -5 GeV -2 . The Cabibbo-Kobayashi Maskawa 
matrix elements: |V CS | « 0.986 and \V c d\ ~ 0.225. 
The masses of heavy mesons: m D + = 1968.30 MeV, 
rn D + = 1869.61 MeV, and m D o = 1864.84 MeV. Isospin 
symmetric masses are used for the light mesons: m v = 
138.04 MeV, m K = 495.67 MeV, and m v = 547.85 MeV 
for the pseudoscalar mesons, and m p = 775.19 MeV, 
Tow = 782.65 MeV, mje = 893.74 MeV, and = 
1019.46 MeV for the vector mesons. The masses of the 
leptons: m e = 0.511 MeV, m /t = 105.66 MeV, and 
m Ve = m V)i = 0 MeV. 
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